Real-time feedback control of a mesoscopic superposition 
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We show that continuous real-time feedback can be used to track, control, and protect a meso- 
scopic superposition of two spatially separated wave-packets. The feedback protocol is enabled by an 
approximate state-estimator, and requires two continuous measurements, performed simultaneously. 
For nanomechanical and superconducting resonators, both measurements can be implemented by 
coupling the resonators to superconducting qubits. 
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Methods for continuous-time tracking and control of 
continuous-variable quantum systems have to-date been 
restricted to localization about single points in energy 
or phase space, a direct quantum analogue of classical 
feedback control [IHl] . One of the most striking features 
of quantum mechanics is the ability for systems to ex- 
ist in spatially separated superpositions, and an impor- 
tant question is whether there is a potentially feasible 
way to use feedback to track and control such highly 
non-classical states, while preserving and protecting the 
coherence. The method for monitoring the square of po- 
sition of a nano-mechanical resonator presented in [B], 
along with the observation in [7 that this measurement 
leaves evolving mesoscopic superpositions intact, have 
provided one of the tools necessary for realizing this task. 

For feedback to be able to preserve a superposition, it is 
not enough for the measurement to merely avoid destroy- 
ing the coherence, it must also replace information that 
environmental noise removes. Although we will show 
that measurements of are not sufficient to do this, 
we can complement them with a quasi-continuous par- 
ity measurement, using a readily available off-resonant 
interaction, and this provides the missing information. 
The central element of this parity measurement was in- 
troduced in |8| in a cavity-QED setting: in our solid-state 
setting it is not only more practical, but we are able to 
implement it in a quasi-continuous, rather than a dis- 
crete manner. The second element essential to feedback 
control is an approximate state-estimator that can track 
the superposition in real-time. Using the estimator we 
develop below, we show that the state can be tracked, 
and the separation of the superposition controlled using 
the feedback cooling method introduced in [51 |S] . 

We note that a previous method has been devised for 
protecting coherence in optical cavities. In a sequence of 
papers the group of Tombesi refined a feedback procedure 
in which a control system is correlated with the parity of 
the cavity field, and if a photon has been lost it is replaced 
by an adiabatic transfer from a single atom [Sl fTUHT^ (see 
also [13 )• This procedure is a very elegant example of 
coherent feedback control [T4UI6] . While their protocol 
could be a useful tool in mesoscopic systems, it is quite 



distinct from our protocol here; it does not provide the 
observer with detailed phase-space information (and is 
thus restricted to preserving only the parity), or provide 
an efficient method to track the state. We note also that 
our protocol does not employ a photon replacement op- 
eration. This is partly because our method separates 
the task of preserving coherence from that of manipulat- 
ing the parity of the superposition (symmetric or anti- 
symmetric): the protocol performs the former, but the 
not the latter. The protection of coherence is our goal 
here since this is the key quantum feature of the state. 

The system that we wish to control is a single quan- 
tum oscillator (with frequency w, position x, and momen- 
tum p) in a superposition of two wave-packets that are 
well-separated in phase space. We will focus on nano- 
mechanical resonators, since preparing these in meso- 
scopic superpositions is a present goal of experimental 
work [T7l421j . It was shown in [7] that a mesoscopic su- 
perposition will remain preserved by a measurement of 
x^, so long as the position-space probability density for 
the superposition is symmetric in x. Further, a contin- 
uous measurement of x^ will create such symmetric su- 
perposition states from an initial thermal mixture. This 
latter fact implies that an x'^ measurement must be able 
to extract parity information (at least for small phonon- 
number) in order to purify the initial mixture into a well- 
defined superposition. 

We begin by deriving an approximate state-estimator 
for the x'^ measurement, which will also reveal how and 
when this measurement will extract parity information. 
This is achieved by constructing an ansatz for the state of 
the resonator, general enough to provide a good approxi- 
mation throughout the evolution. We choose this ansatz 
to be a mixture of symmetric and antisymmetric super- 
positions of pure Gaussian wave-packets. If we denote by 
\G±) the pair of general pure Gaussian states with cen- 
troids {±x,±p), variances Vp and Vx, and symmetrized 
covariance C, then the superpositions 

|±) = (|G+)±|G_))/v/2(T±^ (1) 

are respectively a symmetric state (containing only an 
even number of photons) and an anti-symmetric state 
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(containing an odd number). Here x = R-e[(G+|G_)] is 
the overlap between the two Gaussians. Our full ansatz 
for the state of the resonator is then 

P = ^l+>(+l + (l-^)|-)(-|, (2) 

and contains the six parameters x,p,Vx,Vp,C and P. 
This ansatz is based on the fact that the source of de- 
coherence (the thermal bath) can be described in terms 
of phonon emissions and absorptions [lOj . These change 
the state by one photon (toggle the symmetry), and cor- 
respond to operations linear in x and p, so individually 
they preserve the ansatz. But since the observer does not 
have access to the times at which these random phonon 
events occur, he or she must average over all possible se- 
quences of them. This averaging procedure no longer 
preserves the ansatz, because different sequences shift 
the wave-packets by different amounts, resulting in more 
terms in the density matrix. Nevertheless, so long as the 
measurement is strong enough compared to the rate of 
decoherence (in any case required for effective control), 
the even and odd components will remain nearly pure, 
and so the ansatz should provide a good approximation 
in this regime of "good control" . 

The stochastic master equation (SME) that describes 
the continuous measurement of x'^ is [22J 

dp = -k[M, [M,p]]dt + V2k{Mp + pM-2{M)p)dV, (3) 

with M = x^. The stochastic increment dV is obtained 
by the observer from the continuous stream of measure- 
ment r esults, r{t) = 4:k{M)dt + V2kdW, via the relation 

dV ^ dr/V2k - V8k{M)^dt. (4) 

If the observer is solving the exact master equation, then 
her estimated value of (M) is simply (M)o = Tr[pAf], 
and dV = dW. If the observer is solving an approxi- 
mate estimator, then (M)o is the corresponding approx- 
imate estimate. We are assuming here that our continu- 
ous measurements are quantum- limited, and this requires 
quantum-limited amplifiers. While such amplifiers are 
not yet available for mesoscopic systems, recent progress 
in this area has been rapid, and it appears likely that 
such amplifiers will be realized in the near future |23j . 

To derive the approximate state-estimator, an equa- 
tion of motion for the six parameters, we substitute the 
ansatz into the SME, Eq.(|3|, and assume that the evolu- 
tion preserves the Gaussian form of \G±). The resulting 
approximate state-estimator is 

dx = {p/m)dt + ^xVa^dV (5) 

dp ^ —muj^xdt + ^xCdV (6) 

dV^r = [2C/m ~ i^x?Vl\ dt + ^V^dV (7) 

dC = [Vp/m - muj^V^ - C^x^V^C] dt + ^V^CdV (8) 
dVp = -2mu}^Cdt + ([C^ - tf/i] dV 

(T4 + 2x^)- C^x^] dt (9) 

dP = Vsk {{x% - P{1 - P)dV, (10) 



where we have defined ^ = 2-\/8fc, and the quantities 
(x^)j_ are the estimated expectation values of x^ for the 
states |±) . So long as the wave-packets are well-separated 
(X = 0) we have {x^)l = {x'^)'L =V^+ x? . 

Note that it is the mixing probability, P, that deter- 
mines the level of coherence: P = 1 and P = describe 
perfect coherence, and P = 1/2 denotes zero coherence. 
To preserve coherence the feedback protocol must keep 
the absolute value of the parity, \P\ = |2P — 1|, close to 
unity. Examining the equation for P, we see that this is 
exactly the equation for a classical continuous measure- 
ment that distinguishes between the two states of a single 
bit. The rate at which this information is extracted can 
be read off directly, and is ((a;^)+ - (x^)!) . This 
tells us that the measurement provides parity informa- 
tion so long as (x^)^ ^ {x"^)-- This is true only if the 
Gaussian states in the superposition are close enough in 
phase space to significantly overlap. Our position mea- 
surement will effectively track the location of the wave- 
packets, but will not purify a mixture of the even and 
odd superpositions once they are well-separated in phase 
space. In fact we conclude that no practicable contin- 
uous symmetric measurement of position can determine 
the parity of these states — while the packets cross each 
other (in position space) as they pass x = every half 
period, the spatial period of the oscillations in the over- 
lap become increasingly rapid, so that impractically fine 
resolution in position would be required to distinguish 
odd from even. 

The thermal bath destroys parity information, reduc- 
ing \V\ and thus the coherence over time. The mea- 
surement must continually replace this information to 
preserve the coherence. To solve this problem we now 
show how to make a quasi-continuous parity measure- 
ment, which can be performed in conjunction with the 

measurement. To do this we couple the resonator off- 
resonantly to a single superconducting qubit I24| . The 
interaction is given by i?int — hgdzO^a, where az is the 
Pauli operator for the qubit. Let us denote the eigen- 
states of the ax operator as |-|-) and |— ). If we start the 
qubit in state |+), and switch on the interaction for a 
time r = 7r/g, then the qubit is flipped to state |— ) if and 
only if the resonator has an odd number of photons. This 
correlates the qubit perfectly with the resonator parity. 
If we now make a continuous measurement of <Jx with 
strength p, — which is described by Eq.([3| with M = Ux 
and k ~ p — then it translates to a measurement of the 
resonator parity while the correlation is maintained. The 
thermal damping of the resonator slowly de-correlates the 
two systems, and we must take this into account when 
determining the evolution of P. If we denote the prob- 
ability that the qubit is in state |-|-) by P+, then the 
continuous measurement of the qubit reduces to 

dP+^^P+{l-P+)dU, (11) 

where dU is obtained from the measurement record. 
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Txit), by dU = drxit)/^/2iJ.— ^/8fi{ax)dt. In our case, the 
degree of correlation between the qubit and the parity is 
fully determined by the conditional probability for the 
resonator to have positive parity, given that the qubit 
is in state |+). Denoting this conditional probability 
by 1 — /?, the observer's estimated mixing probability is 
P = (1 - i3)P+ + /3(1 - P+). Since the thermal damp- 
ing de-correlates the qubit, we reset and re-correlate the 
qubit at regular intervals. We obtain good tracking of the 
parity if we allow the estimator to assume that during the 
correlation step, P decays due to the thermal damping, 
but at the end of the step, (3 = (perfect correlation), so 
that the estimator can reset P+ to the current value of P. 
Between re-correlations, the evolution of /3 is determined 
by the thermal damping, and we give this below. We 
do not take into account relaxation (or dephasing) of the 
qubit in the above analysis, because simple quantum (or 
in fact classical) error correction can be used on the qubit 
to eliminate the effects of relaxation on the measurement. 

Both the and measurements must be switched- 
off while the qubit is correlated with the resonator parity, 
since they will interfere with it. We thus make the cor- 
relation step short compared to the effective frequency 
of the resonator, which is feasible since the effective fre- 
quency can be relatively low (e.g. 250 kHz, see below). 

Before we simulate our state-estimation procedure, we 
need to know how the thermal damping of the resonator 
changes our six variables, so as to include this dynam- 
ics in the state-estimator. These equations of motion 
are necessarily an approximation, because as mentioned 
above the damping does not preserve the ansatz. In 
the well-separated regime (x = 0), our approximate 
damping equations are x = —(7/2)0;, p = —(7/2)^, 
Vx = -liVx-Vj"-), Vp = -7(T/p-y7'^), and C = -2jC. 
In these equations 7 is the damping rate of the resonator. 



and Vj^ = {2nT + 1) and V^^ = {2nT + 1) are the 
variances of the harmonic oscillator at temperature T. 
The quantities A^ = h/{2muj) and A^ = hmuj/2, are, re- 
spectively, the variances of x and p in the ground state. 
The parameter ut is the average number of phonons that 
the resonator would have if it were at the bath tempera- 
ture T. The evolution of the qubit /resonator correlation 
between correlation steps is 



/3 



7 Vx + x + Vp+p-1 {(3-1/2), 



(12) 



where the tildes indicate that the means and variances 
are those of the dimensionless variables x = {a + a^)/\/2 
and p = —i{a — a')/y/2. Note that the expression in the 
square brackets is merely the estimated value of 2{a^a). 

We now perform simulations to verify that our pro- 
tocol is able to track and preserve the superposition, 
as well as control the phase-space separation of its two 
components. We do this using experimentally realistic 
parameters, as follows. Nano-resonators typically have 
frequencies / = lli/{2tt) in the range 10 — 100 Mhz. 
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FIG. 1. (color online) Time series' for the true and estimated 
values of the evolution of various observables of the resonator 
under the feedback control protocol, (a), (b) and (c) have 
resonator quality factor Q = 200, and temperature T = 0. (d) 
and (e) have Q = 500, and T = 20 mK (nr = 12). The dark 
(blue) solid lines are the true values, and the Ughter soUd lines 
(red or blue) are the estimated values. In (c) the dashed line 
indicates when the feedback is heating (high value), dormant, 
or cooling (low value) the system. 



By modulating at frequency A the coupling of the res- 
onator with the probe qubit mediating the x^ measure- 
ment, we can reduce the oscillation frequency (from the 
point of view of the measurement and thus the observer) 
to V = a; — A. This means that the state-estimation 
need only work on the timescale of rather than w. 
We chose uj/{2tt) = 50 MHz, and v/{2tt) = 250 kHz. 
Since nano-resonators typically have Q factors of 10^, 
the damping rate is 7/(27r) = 5 x 10'^ s^^, and this gives 
an effective Q of 1^/7 = 500. Both the x^ measurement 
and the energy interaction Hi^t are obtained by using an 
off-resonant (perturbative, 2"''-order) interaction with a 
qubit [7j . The maximum rate of direct coupling between 
qubits and resonators is ^ 10* s~^ [5S], so that of 2"''- 
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order interactions is 



10^ 



To fix the correlation 



step at 1/20 of the effective resonator period, we need 
the interaction rate g = \2.^v sa 8 x 10^ s"^ . The rates 
of the and Uz measurements should ideally be much 
larger than ^tit, and the former should be less than or 
v to avoid undue heating. For the following simulations 
we chose k = v /[2QQ-k) and two values for ^, = v and 
/i = v/2. 

To simulate the continuous measurements of and cr^ 
in the presence of thermal damping we use the quantum 
Monte-Carlo method described in We also use the 
approximate state-estimator to implement the feedback 
cooling method devised in [21 [9] to stabilize the amplitude 
(separation) of the superposition. The cooling method 
involves raising the potential (increasing the frequency 
of the oscillator) when the estimated value of (a;^) is at 
a minimum, and reducing the frequency when it is at a 
maximum. To determine when (a;^) has reached an ex- 
tremal point, we fit a quadratic to a smoothed version of 
the time-series for the estimated (a;^), looking back for 
one twentieth of a period. In our simulations we stabi- 
lize the state by applying cooling (heating) — the lat- 
ter using the reverse algorithm — when the estimated 
average phonon number, (ri)c, is greater than 25 (less 
than 16). This feedback stabilization is turned off when 
16 < (n)o < 20. 

In Fig. [T](a)-(c) we simulate the estimation and con- 
trol process using the parameters above, but with a 
lower quality resonator {Q = 200), and with the tem- 
perature set to zero. The initial state is an equal mix- 
ture of the even and odd superposition states, with 
(x,p,14, Vp,C) = (6V2A^,0, A2, A2,0). Here we start 
the estimator with the correct values of the parameters, 
but we have verified that if the initial guess is wrong, 
the estimator locks on in two or three oscillation pe- 
riods. The estimator appears to be quite robust, but 
to increase robustness our estimator truncates the mea- 
surement noise if it is outside 4 standard deviations, and 
resets the variances if one of the following conditions oc- 
curs: 14 < 0,yp < Q.V^Vp < 0.9(C2 -f 0.25fi2),T4 > 
12A2, or Vp > 12A2. In Figs.[l];a), (b), and (c) we show 
the true and estimated values of, respectively, the ab- 
solute value of the scaled position, parity, and average 
phonon number as functions of time, for a single run. 
These plots show that the estimator tracks the location 
and parity of the state well, and shows the degree to 
which the energy is stabilized by the feedback. The par- 
ity time-series also shows that the measurement is ef- 
fective at keeping the state highly pure, and exhibits the 
quantum jumps in parity. These jumps are due to the in- 
terplay of the strong quasi-continuous measurement and 
the much weaker thermal decoherence (phonon damping) 
— the measurement is able to reveal the phonon emission 
events that flip the parity. 

In Fig.[l];d) and (e) we simulate the process again, but 
this time with Q = 500, and set the temperature to be 



20 mK. For our 50 MHz resonator this means that the 
thermal occupation number is tit = 12. We see that the 
estimator still tracks well the position and the parity. 
But this time, due to the increase in ht, the measure- 
ment is not able to maintain the purity as effectively as 
before. Nevertheless, the estimator still reveals the quan- 
tum jumps in the parity due to the thermal environment. 
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